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Abstract 

The space- time conservation element and solution 
clement (CE/SE) method is used to study the sound- 
shock interaction problem. The order of accuracy of 
numerical schemes is investigated. The linear model 
problem governed by the 1-D scalar convection equa- 
tion, sound-shock interaction problem governed by 
the 1-D Euler equations, and the 1-D shock-tube 
problem which involves moving shock waves and 
contact surfaces arc solved to investigate the order 
of accuracy of numerical schemes. It is concluded 
that the accuracy of the CE/SE numerical scheme 
with designed 2nd-order accuracy becomes 1st or- 
der when a moving shock wave exists. However, the 
absolute error in the CE/SE solution downstream 
of the shock wave is on the same order as that ob- 
tained using a fourth-order accurate essentially non- 
oscillatory(ENO) scheme. No special techniques arc 
used for either high- frequency low-amplitude waves 
or shock waves. 

1. Introduction 

In the computational acroacoustic problem, the 
propagation of disturbances of small amplitude 
needs to be captured. Thus a stringent requirement 
is placed on a numerical algorithm. And further, 
when a shock wave and acoustic waves exist at the 
same time, it becomes more challenging for a nu- 
merical simulation. In the current field of CAA, 
the high-order accuracy finitc-diffcrcncc method is 
popularly used for the propagation of high-frequency 
low-amplitude waves because a traditional 2nd-ordcr 
accurate scheme is not accurate enough for this kind 
of problem. When a shock wave exists, the numer- 
ical scheme also has to be able to accurately cap- 
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turc the shock wave without oscillations. Tradi- 
tional high-order accurate shock-capturing methods 
are classified into linear and nonlinear methods. The 
EXO schemes belong to the latter class. The order 
of accuracy of an ENO scheme is investigated for a 
sound-shock interaction problem in [1]. A very in- 
teresting conclusion drawn in [1] is that the accuracy 
of the solution downstream of the shock wave drops 
to first order if the shock wave is not located at a 
mesh point. The designed accuracy downstream of 
the shock wave can be achieved by using subcell res- 
olution in a numerical algorithm such that the shock 
wave is located exactly on a mesh point. However, 
such a strategy is not practical for general multi- 
dimensional problems, since the shock wave could be 
curved and moving in the entire domain. Therefore, 
the advantage of using high-order accurate methods 
in the study of unsteady flows with shocks is ques- 
tionable. 

The space-time CE/SE method is an innovative 
numerical method for solving conservation laws. It 
is different in both concept and methodology from 
the well-established traditional methods such as the 
finite difference, finite volume, finite element and 
spectral methods. It is designed from a physicist’s 
perspective to overcome several key limitations of 
the traditional numerical methods. 

Simplicity, generality and accuracy are pursued in 
the development of the CE/SE method. Its salient 
properties arc summarized briefly as follows. First, 
the concepts of conservation element and solution 
clement arc introduced to enforce both local and 
global flux conservations in space and time instead 
of in space only. Second, all the dependent variables 
and their spatial derivatives are considered as indi- 
vidual unknowns to bo solved for simultaneously at 
each grid point. Third, no approximation techniques 
other than Taylor’s scries expansion, no monotonic- 
ity constraints, and no characteristic-based tech- 
niques arc used in the design of the scheme. A de- 
tailed description of this method and the accompa- 
nying analysis arc set forth in [2-4]. 
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A variety of numerical tests have been per- 
formed previously to illustrate the robustness of this 
method. For the CE/SE Euler solver, highly ac- 
curate numerical solutions have been obtained for 
various flow problems involving discontinuities, such 
as shock waves, contact surfaces and even their in- 
teractions [5.6]. Moreover, applications of the same 
Euler solver to computational acroacoustics (CAA) 
problems reveal that the accuracy of the results 
is comparable to that of a 4th-ordcr compact dif- 
ference scheme even though the current solver is 
only 2nd-ordcr accurate. Further, the non-reflecting 
boundary conditions can be implemented in a simple 
way without involving characteristic variables. The 
solver can be applied to subsonic, transonic, and su- 
personic flows in the same form without using the 
characteristic-based techniques. Results show that 
the present solver can handle both continuous and 
discontinuous flows very well [7-14]. 

In this paper, the test problems used in [1] are 
solved by using the CE/SE method. Problems gov- 
erned by the linear convection equation and the Eu- 
ler equations arc studied. The order of accuracy of 
the CE/SE numerical schemes is investigated. Nu- 
merical solutions are compared with those obtained 
by using the ENO scheme. 


2. Numerical Test Problems 


2.1. A linear model problem. 
Consider the scalar equation 


du du n 

a + "aF = ° 


(i) 


where the wave speed a is 


a = 


{ 


2 . 

E 


x < x s 
X > x s 


(2) 


where x s = 0.5 given in [1]. The initial conditions 
are described as 


u{x> 0) 


{ 1/2. x < x 8 , 
l, x > x s 


( 3 ) 


The domain is 0 < x < 1. The inflow boundary 
condition is 


e — 0.5, a = 1. and Courant number 0.8. The de- 
tails of the scheme arc given in [2] and won’t be 
repeated here. 

First, an even number cells with 64, 128. 256. 512. 
1024 uniform mesh intervals is used. The computed 
perturbation u f (x, t) = u(*r. 0 - u(x.0) at t — 107a 
and its logio error for different cell numbers arc 
shown in Fig. 1. It can be seen that the scheme 
is 2nd-ordcr accurate in the entire domain except at 
the shock location. In Fig. 2, the same results arc 
shown for an odd number cells with 65.129,257.513 
uniform mesh intervals. It is still 2nd-ordcr accurate 
in the entire domain since the space-time CE/SE 
method uses a staggered mesh. There is a grid point 
located at x 6 = 0.5 no matter whether an even or 
an odd number cells is used. An EXO-4-3 scheme, 
which is fourth-order accurate in space and third- 
order accurate in time, is used in [1]. For this linear 
problem, the accuracy of the ENO-4-3 solutions re- 
mains 4th order for an even number cells because 
there is a mesh point at the shock location, while 
the accuracy drops to first-order downstream of the 
shock wave for an odd number cells because the 
shock wave is within a cell. 

Second. x 6 = 0.5 is replaced by x s = 0.5 + ds with 
ds — 0.00035, which results in no mesh points co- 
inciding with the shock location even on the finest 
mesh used. The solution is shown in Fig. 3. The ac- 
curacy downstream of the shock wave drops to first 
or zeroth order. It was found that the dominant 
error is from the phase error due to ds which can 
not be resolved. From this test result, it can be 
concluded that the designed 2nd-ordcr accuracy can 
not be achieved if a shock is within a cell for the 
CE/SE scheme applied to linear problems. In order 
to obtain more relevant conclusion about the perfor- 
mance of the CE/SE method, the physical problem 
governed by the Euler equations is considered next. 


2.2. Shock- Sound Interaction 

The shock-sound interaction problem described in 
[1] is governed by the 1-D Euler equations, whose 
conservative form is expressed as 


du m df m 

dt dx 


777 = 1. 2, 3 


(5) 


in which 


u(Q.t) = (1 — £s\mjt)/2 (4) 

where s = 0.001 and u; — 8tt. The period of the 
acoustic wave T\ is define as 2 tt/u;. The outflow 
boundary condition is imposed by using the ana- 
lytical solution given in [1]. The ID CE/SE a-e-a 
scheme constructed for Eq. (1) is used here with 


u i = P, u 2 = p\\ u 3 = p/(- - 1) + pu 2 / 2 (6) 

and 

h = U2 (7) 

h = (7 - 1)«3 + (3 - -)(t*2) 2 /(2«i) (8) 

h = - u 2 u 3 /wi - (l/2)(7 - 1 )M 3 /(«i) 2 (9) 


2 

American Institute of Aeronautics and Astronautics 


where p,v,p. and 7 are the mass density, velocity, 
static pressure, and constant specific heat ratio, re- 
spectively. 

The spatial domain is 0 < x < 1. The main flow 
is from left to right. A shock wave is initially lo- 
cated at x = 0.5. The initial conditions of main 
flow variables at the left side(upstrcam) and right 
side (downstream) of the shock wave arc described 
as follows: 

(p. i\p) L = (1.0. 2 . 0 . 1 . 0 / 1.4) x < 0.5 ( 10 ) 

(p. v.p)n = (2.6666. 0.75. 4.5/1. 4) x > 0.5 (11) 

The acoustic disturbance is introduced at x — 0. 
The flow variables at the inlet are defined as 

p( 0 ,() =pl(1 -esinwf) ( 12 ) 

P(0,t) = pl\p( 0. t)/p L ) lh (13) 

i'(0. t) = v L - (<r(0- 1)- cl) (14) 

7-1 

where u) = 36-, e = 0.001, and c = y/jp/P being the 
local sound speed, thus cj_ = \J~ Pl/pl is the sound 
speed upstream of the shock wave. At the outlct(x — 
1 ). the non-reflecting boundary condition is imposed 
to let flow propagate out the computational domain. 

In this problem, the shock wave moves around the 
initial location due to the interaction with the acous- 
tic wave. Thus, the shock wave is not located at 
a mesh point. Therefore, the first-order accuracy 
was achieved downstream of the shock wave using 
the EXO-4-3 scheme, which is shown in Fig. 4. 
However, in this problem, the shock wave docs not 
move outside the cell in which it is initially located. 
Thus, inclusion of subcell resolution in the ENO- 
4-3 scheme can be done within this cell to exactly 
resolve the shock location. The 4th order accuracy 
downstream of the shock wave was achieved, which is 
shown in Fig. 5. Subcell resolution can be achieved 
more easily for the linear model problem described 
in 2 . 1 , since the shock wave does not move and re- 
mains at x = 0.5. However, if shock waves arc mov- 
ing across the cells as the solution evolves, the use of 
subccll resolution in a scheme would be impractical 
due to its cost and complexity. 

This problem is solved here again using the ID 
CE/SE Euler solver. The same initial conditions 
and the boundary condition at the inlet are used. 
At the outlet (x = 1 ). the non-reflecting boundary 
condition which is not based on the characteristic 
theory is imposed as 

(Um)] = (U m (U mx )] = 0 (15) 


The boundary condition at the outlet is not men- 
tioned for EXO scheme in [ 1 ]. Two parameters e 
and a for controlling numerical dissipations in the 
CE/SE Euler solver arc set as 0.5 and 1. and the 
Courant number is 0.8 in the calculations. The com- 
puted acoustic wave solutions p f (x,t) = p(x.t) - 
p(x. 0) and p f {xJ) = — p(x, 0) at t = 307\ arc 

shown in Fig. 6 along with the exact solutions. The 
logio errors of both //(x, t) and p'(x. £) arc shown in 
Fig. 7. It can be concluded that the accuracy of the 
CE/SE solution is first order in the entire domain 
except at the shock location. The absolute error of 
the CE/SE solution downstream of the shock wave 
is on the same order as that of the EXO-4-3 solution 
shown in Fig. 4 using the same 512 uniform mesh 
intervals. 

2.3. 1-D Shock- Tube Problem 

Finally, the 1-D shock-tube problem named Socks 
problem [ 2 ] is solved here in an accuracy study us- 
ing the same solver mentioned in 2.2. In this prob- 
lem, the shock wave and contact surface propagate 
through the entire spatial domain during the time 
interval of computation. The numerically computed 
solutions of pressure and density at t = 0.24 with 
CFL — 0.44 arc shown in Fig. 8 along with the ex- 
act solution. The logio errors of pressure and density 
arc shown in Fig. 9. respectively. It can be seen that 
the solution is first-order accurate in the entire do- 
main except at the shock wave and contact surface 
locations. No solution obtained by an ENO scheme 
is available in [ 1 ], 

3. Conclusions 

Three test problems have been solved to inves- 
tigate the accuracy of the CE/SE method for un- 
steady compressible flows with shock waves. Gener- 
ally speaking, first-order accuracy can be obtained 
for the CE/SE method without using any special 
techniques for either high-frequency low-amplitude 
waves or shock waves. The absolute error in the 
CE/SE solution downstream of the shock wave is 
on the same order as that obtained by using fourth- 
order accurate ENO scheme for sound-shock inter- 
action problem. It can be concluded that CE/SE 
method can produce accurate solutions for acroa- 
coustic problems involving shock waves in a simple 
way. 
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Figure 1: The CE/SE solution at t = 107a (upper) 
and the point wise error distribution using even num- 
bered cells for a discontinuity located at x — 0.5. 


Figure 2: The CE/SE solution at t — 107a (upper) 
and the pointwisc error distribution using odd num- 
bered cells for a discontinuity located at x = 0.5. 
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Figure 3: The CE/SE solution at t = 107\ (up- 
per) and the point wise error distribution using even 
numbered cells for a discontinuity located at x = 
(bo -t- 0.00035. 


log ! Error ( 



Figure 4: The ENO-4-3 pointwisc error of pressure 
for the sound-shock interaction problem at t — 30T\. 


log! Error I 



Figure 5: The EXO-4-3-SR pointwisc error of pres- 
sure for the sound-shock interaction problem at t = 
30 T x . 
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Figure 6: The CE/SE solutions of the sound-shock 
interaction problem at t = 307V 




Figure 7: The CE/SE pointwisc errors for the sound 
shock interaction problem. 
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Figure 8: The CE/SE solutions of the shock-tube Figure 9: The CE/SE point wise errors for the shock- 
problem (Sod’s problem) at t — 0.24 with CFL — . tube problem with CFL = 0.44. 

0.44. 
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